Quantum transport in the presence of a finite-range time-modulated potential 
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Quantum transport in a narrow constriction, and in the presence of a finite-range time-modulated 
potential, is studied. The potential is taken the form V(x, i) — Vo 6(x)8(a — x) cos(cji), with a the 
range of the potential and x the transmission direction. As the chemical potential fi is increasing, 
the dc conductance G is found to exhibit dip, or peak, structures when fi is at nhto above the 
threshold energy of a subband. These structures in G are found in both the small a (a <C Af) and 
the large a (a S> Xf) regime. The dips, which are associated with the formation of quasi-bound 
states, are narrower for smaller a, and for smaller Vo- The locations of these dips are essentially 
fixed, with small shifts only in the case of large Vo- Our results can be reduced to the limiting case 
of a delta-profile oscillating potential when both a and Voa are small. The assumed form of the 
time-modulated potential is expected to be realized in a gate-induced potential configuration. 
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PACS numbers: 72.10.-d, 72.40.- 



I. INTRODUCTION 



Inelastic scattering processes in quantum transport 
have drawn continuous attentions in the recent past. 
One of the common model used is to invoke a time- 
modulated potential, with a certain spatial profile, to 
the systemi 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 This model was first utilized 
by Biittiker and Landauer— who, in considering a time- 
modulated barrier, attempted to extract a traversal time 
of a particle tunneling across the barrier. Later work 
from other groups have considered inelastic scatterings 
in different configurations such as the double oscillat- 
ing barriers 3 ", the time-modulated barrier with a delta 
profile 4 ^, and the oscillating quantum well in between 
two static barriers 7 ^. The model has also been ex- 
tended to incorporate the inelastic effects due to phonons 
by introducing a time-modulated potential involving the 
phonon operatorsiii* 1 ^* 1 ^* 1 ^ In all these studies, the fea- 
tures arising from the inter-mixing between the elastic 
scattering, through a static potential, and the inelastic 
scattering, through a time-varying potential, is empha- 
sised. In particular, these studies have demonstrated, 
among others, the interesting feedback effect of the in- 
elastic scattering on the elastic channel. Even though the 
above model is appropriate only for inelastic processes 
that preserve the phase coherence of the transmitting 
particles, the model has practical importance because 
the coherent inelastic scattering (CIS) can be realized, 
at least, in the case when the time- modulated potential 
is well-specified. 

A possible realization of CIS processes in nanostruc- 
tures is expected to be found in gate-controlled quan- 
tum point contacts (QPC), as shown in Fig. ^ Similar 
gate-induced potential configuration has been suggested 
by Gorelik et al.^ who considered microwave- induced 
effects on the Josephson current through a narrow con- 
striction (NC). Their focus is on the resonance of the 
microwave frequency with the energy levels of the An- 
dreev bound states formed in the NC, which has both 
ends connected to superconducting electrodes. For our 



purposes here, a simple exhibition of the CIS is expected 
to be found readily in a normal state gate-controlled 
QPC. Recent development in the split-gate technology 
has made possible the fabrication of such gate-controlled 
QPCi 15 ' 16 The split-gates, when negative biased, define 
electrostatically a NC on a two-dimensional electron gas 
(2DEG). The dc quantum transport properties of these 
QPC systems has been studied intensively? 1 ^*^ More re- 
cently, there are growing interest in the time-dependent 
properties, such as the effects of photons, in these QPC 
systems A^^SSiSiiSSiS It is thus legitimate to consider the 
quantum transport in a NC which is acted upon by an 
additional, and ac biased, gate, as shown in Fig. 2] This 
ac biased gate, which is different from the split-gates that 
define the NC, induces on the NC a time-modulated po- 
tential. The scattering of the conduction electrons by this 
time-modulated potential is both coherent and inelastic. 

There is another reason why quantum transport in 
a NC in the presence of a time-modulated, and gate- 
induced, potential is interesting. This is closely re- 
lated to the density of states (DOS) structures in the 
NC. The energy levels in the NC are quantized into 
one-dimensional subbands so that the DOS is singular 
at the subband bottoms. In the presence of attrac- 
tive impurities, such singularities in the DOS lead to 
dip structures in the dc conductance G, as the chem- 
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FIG. 1: Sketch of the gated QPC in which a NC is connected 
adiabatically at each end to a 2DEG electrode. The gate 
induced a finite-range time-modulated potential in the NC. 
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ical potential p is increasing, 24 i 25 . 26 i 27 i 28 i 29 . 30 i 31 i 32 The 
dip structures occur when u is just below a subband 
edge. According to Bagwell^ these dip structures are 
associated with the formation of impurity-induced quasi- 
bound-states (QBS)i25( The wavefunction at this energy /i 
and in this subband is evanescent along the transmission 
direction. Hence for the case of an attractive impurity, 
a QBS splits off from each subbandi^ An electron orig- 
inally in a propagating state in other subband can thus 
be scattered elastically into and be trapped by this QBS. 
This gives rise to dip structures in G. 

The QBS features in G are found also when a point 
barrier oscillates in a NC£ In this case, for a not-too- 
large oscillation amplitude, the dc conductance G ex- 
hibits dip, or peak structures when fi is at nhuj above 
a subband edge. These structures correspond to the 
situation when the electrons can make transitions, via 
inelastic processes, to the QBS just below the subband 
edge. That there is QBS, induces by the point oscillat- 
ing barrier, below each subband edge is demonstrated 
by Bagwell^ from the energy poles in the current trans- 
mission coefficients. The imaginary part of the energy 
poles are negative, which is consistent with the nature 
of the QBS. The existence of these QBS is again due to 
the singular DOS near a subband edge. It is important 
to ask whether such QBS features in G persist in the 
case of a finite-range time-modulated (FRTM) potential. 
This question has not been addressed before and, if the 
QBS features did persist in a FRTM potential, will have 
important implications to time-dependent properties of 
QPC systems. Furthermore, since the potential is ex- 
pected to be gate-induced, the problem is within reach 
of the recent experimental capability. 

In this paper, we simplify the problem by assuming 
the the gate-induced potential in the NC to be repre- 
sented by the form Vo9(x)9(a — x) cos(wt), where a is the 
range of the potential and x is the transmission direc- 
tion. Our simplification is in replacing the smooth longi- 
tudinal potential profile, of which the potential builts up 
within a longitudinal distance of order Ap, by an abrupt 
profile. The abruptness of the profile is expected to do 
nothing except introducing additional multiple scatter- 
ings between the two abrupt edges of the potential. This 
results in additional harmonics in G. Thus for the cases 
when the magnitude of these harmonics is small, our re- 
sults are expected to resemble qualitatively the features 
in a smooth profile potential. An explicit smooth-profile 
consideration, however, is left to the further study. 

Using this small-harmonic-magnitude criteria, we find 
the QBS features in G in both the small a (a <C Af) 
and the large a (a 3> Af) regime. Our results can be 
reduced to the limiting case of a delta-profile oscillating 



potential when both a and Vqo are small. In our calcula- 
tion, the inelastic scattering is solved nonperturbatively. 
From our results, we note that even within the small- 
harmonic-magnitude criteria, the inelastic scattering has 
to be treated beyond one sideband approximation. The 
sideband index n labels those electrons which net energy 
change is nTioj. 

In Sec. II we present the formulation for the inelastic 
scattering and the connection of the current transmis- 
sion coefficients with the conductance G. In Sec. Ill we 
present numerical examples illustrating the QBS features 
in a FRTM potential. Finally, Sec. IV presents a conclu- 
sion. 



II. THEORY 

In this section, the inelastic scattering problem is for- 
mulated and the equations for the current transmission 
and reflection coefficients are obtained. The conductance 
G is then expressed in terms of these coefficients. 

The QPC is modelled by a NC connecting adiabati- 
cally at each end to a 2DEG. Hence the transmission of 
the electrons into, or out of, the NC region is adiabatic. 33 
The gate-induced potential is assumed to affect only the 
NC region of the QPC. Therefore we need only to for- 
mulate the inelastic scattering in the NC region. The 
NC is taken to have a quadratic transverse confinement 
potential lo 2 y 2 . The gate- induced potential is taken the 
FRTM form 

V(x,t) = V 6(x)9(a-x)cos(cjt), (I) 

which connection with a smooth-profile potential has 
been discussed in the previous section. 

Choosing the energy unit E* = h 2 kp/2m*, the length 
unit a* = 1/kp, the time unit t* — h/E*, and Vq in units 
of E* , the dimensionless Schrodinger equation becomes 

[-V 2 + uj 2 y y 2 + V(x, t) ] t) = i- *(x, t). (2) 

Here kp is a typical Fermi wavevector of the reservoir and 
m* is the effective mass. The transverse energy levels are 
quantized, with e n = (2n + I)wj,, and </> n (y) the wave- 
function. The FRTM potential is uniform in the trans- 
verse direction and does not induce inter-subband tran- 
sitions, leaving the subband index n unchanged. Thus 
for a nth subband electron incident along x, and with 
energy fx, the scattering wavefunction can be written in 
the form ^(x, t) = (f> n {y) t), where 3 ^ 
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(3) 



and n, to are the final subband and sideband indices, re- 
spectively. The effective wavevector for an electron with 
energy e and in the nth subband is given by k n (e) = 
\Fe — (2n + l)w y . The sideband index to corresponds to 
the net energy change of mTiuj for the outgoing electrons. 

It is very important to note that had the length of the 
NC been infinite, and the range of the potential V(x, t) 
were extended to cover the entire NC, the longitudinal 
wavevector k n would be a good quantum number so that 
no real transition could have occurred. However, as long 
as V(x, t) has a finite range, k n is no longer conserved, 
and real transitions from k n (e) to k n (e ± mu>) are per- 
mitted for electrons traversing the potential. Thus the 
finiteness in the range of the time-modulated potential 
alone makes possible the absorption of energy by the elec- 
trons for arbitrary u. This picture holds regardless of the 
range, long or short, and the profile, abrupt or smooth, of 
the potential. The mathematical statement of the above 



physical picture turns out naturally, and given by Eq.(4), 
in the following. 

The expressions for the reflection and the transmission 
coefficients can be obtained from matching the wavefunc- 
tions, and their derivatives, at the two ends of the FRTM 
potential. For the above matching to hold in all time, the 
integration variable e in Eq.(3) has to take on discrete 
values fi ± mu>. Hence we can write A n (e) and B n (e) in 
the form 

F n (e) = J2 F n(m)S(e-fi-mu;), (4) 

m 

of which F n (e) refers to either A n (e) or B n (e). After per- 
forming the matching, and the current reflection coeffi- 
cients r n (rn) eliminated, we obtain the equations relating 
A n (m), B n (m), and the current transmission coefficients 
t„(m), 



tn(m) = J2 [A n {m')e- iK -^ m '> + B n (m') e~ iK ^ m ' m >] J m - m ,(V /u), 



(5) 



and 



k n {ii + mu)t n {m) = ^fc n (M + mV) [^(m') e"*- ( m ' m '' 8 - ^(m') e^^" 

m' 

x J m - m '(Vo/uj), 

to, to') + B n (m!) K n (to, to')] J m _ m < (Vq/u), 



(6) 
(7) 



where (to, to') = k n (n+mu))±k n (iJ,+m'uj). Equations 
(5)- (7) can be shown explicitly to reduce to the corre- 
sponding equations for the delta-profile time-modulated 
potential in the a — > limits 

The zero temperature conductance is given by G = 

(2e 2 /^) En=o G n, where G n = £ m G™, and N + 1 is 
the number of propagating subbands in NC for the chem- 
ical potential fi. The contribution to G from electrons 
incident from subband n and transmit into sideband m 
is given by G™, which is related to the transmission co- 



efficients, given by 

G™ = [fc n (/i + to^/M")] I t»(m) | 2 . (8) 

Solving Eqs.(5)-(7), we obtain t n (m), A n (m), and 
B n (m), using which the current reflection coefficient 
r n (m) can be calculated, 

r n (m) = [A n (m') + B n (m')] J m - m /(V /uj) - 6 m0 . 

m' 

(9) 
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We solve the coefficients r„(m), and t n (m) exactly, in 
the numerical sense, by imposing a large enough cutoff 
to the sideband index. The correctness of our procedure 
is checked against the conservation of current condition, 
given by 



k n (p + muS) 
k„ (p.) 



[\t n (m)\ 2 + \r n (m,)\ 2 



1. 



(10) 



III. NUMERICAL RESULTS 

We calculate, in the following, the conductance G of 
a NC acted upon by a FRTM potential. The FRTM 
potential does not induce inter-subband transitions and 
so each occupied subband contributes independently to 
the total conductance. Thus it suffices for our purposes 
here to present the conductance of only one subband, 
which we take it to be the lowest one. 

In this section, the behavior of G with respect to the 
chemical potential fj, is studied. Since G depends also on 
the potential range a, and the oscillating amplitude Vq, 
we present the behavior of G in four situations. Firstly, 
this G behavior is shown for a fixed while varying Vq. Sec- 
ondly, the G behavior is presented for a fixed Vd — Voa, 
while varying a. This is to make connection with the 
delta-profile results, in which the time-modulated poten- 
tial is of the form VdS(x) cos(uit). Thirdly, the G behavior 
for Vo fixed, while a varying, is presented. The fourth sit- 
uation is to compare the G behavior for different u>. In 
addition, we present the time averaged spatial distribu- 
tion for the scattering state, of which the incident energy 
is very close to the QBS structure, and the QBS is about 
to occur. Finally, we compare the one-sideband approx- 
imation results with the nonperturbative results. 

In our numerical examples, the NC is taken to be 
that in a high mobility GaAs-Al^Gai-^As with a typ- 
ical electron density n ~ 2.5 x 10 11 cm -2 , and m* = 
0.067m e . Correspondingly, our choice of energy unit 
E* = h 2 k 2 F /(2m*) = 9meV, length unit a* = l/k F = 
79.6 A, and frequency unit uj* = E*/h = 13.6THz. We 
also take u> y = 0.035, such that the effective NC width 
is of the order of 10 3 A. In the following, in presenting 
the dependence of G on fi, it is more convenient to plot 
G as a function of X instead, where X = [(/j,/ujy) + l]/2. 
The integral value of X is the number of propagating 
channels. 

In Figs. 2 and 3, G is plotted against X for a fixed, 
while Vo is varying. The a is fixed at 150, 62, and 10 in 
Figs. 2, 3(a), and 3(b), respectively. The frequency uj is 
taken to be 0.014, which energy interval u> corresponds 
to an interval AX = uj/(2uj v ) = 0.2 on the ordinate. The 
threshold, or the subband edge, is at X = 1. We note 
that in all the cases shown, a major dip structure occurs 
at X — 1.2, which corresponds to X — AX = 1. This is 
the QBS features because the electron with energy at X 
can make transition to the subband edge by giving up an 
energy u>. We note also that, in general, for larger Vq, the 
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FIG. 2: Conductance G as a function of X for a = 150, and 
ui = 0.014. The potential oscillating amplitudes are: "a" 
Vo = 0.016, "b" Vo = 0.012, "c" V = 0.008. The curves 
are vertically offset for clarity. 



structures at X = 1 + N AX are more evident. These 
are the situations when the electron can emit an energy 
of Nuj and makes transition to the subband edge. The 
wavelength of the electron decreases as X is increasing. 
The relation is given by A = 2u / yj2uj y {X — 1). At the 
location of the first dip, when X = 1.2, we have A = 53. 
Thus Fig. |21 corresponds to the case of a long potential 
range, with a ~ 2.8A, and Fig. 3(b) corresponds to the 
case of a short potential range, with a ~ 0.19A, near the 
occurrence of this first dip. 

Besides the QBS features, there are harmonic struc- 
tures in Fig. [21 and in Fig. Efa). The structures are 
smaller for the lower V~ - That these harmonics are asso- 
ciated with the multiple scattering between the abrupt 
edges of the potential can be identified from a resonance 
relation: A = 2a/n, with n & positive integer. Corre- 
spondingly, the harmonic peaks are at X n = 1 + AX n , 
with AX n = (nir / a) 2 / (2uj v ) . According to the above 
estimate, the first five harmonic peaks are at X n ~ 
1.006, 1.025, 1.056, 1.1, and 1.16, which correspond quite 
reasonably to that in Fig. 2. However, for X > 1.2, 
the harmonic peaks correspond more closely to X = 
1.2+AA„. This can be explained as follows: the harmon- 
ics for A > 1.2 are contributed mostly from those elec- 
trons that give away an energy of uj so that the harmonics 
at X are very similar to that at X — uj, or X — 0.2 for our 
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FIG. 3: G as a function of X, for u = 0.014. (a) a = 62, Vo = 
0.003,0.012,0.024; (b) a = 10, Vo = 0.020,0.028,0.036; (c) 
FRTM potential result for a = 10, Vo = 0.036 (solid curve). 
Delta-profile result for Vd = 0.36 (dashed curve). 
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FIG. 4: G as a function of X for V d = 0.2, and lu = 0.014. 
The parameter a = 10 (solid curve); 20 (dash-dotted curve); 
30 (dotted curve). 
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FIG. 5: G as a function of X for V = 0.02, and w = 0.014. 
The parameter a = 10 (solid curve); 20 (dash-dotted curve); 
30 (dotted curve). 



cases. In Fig. 2, we see that the harmonic amplitudes 
are essentially smaller than the dip structure at X = 1.2. 
Using our small-harmonic- amplitude criteria, the QBS 
features are expected to be evident in a smooth-profile 
FRTM potential. Similar arguments can be applied to 
Fig. 3(a) to establish the harmonic peak locations, but 
we do not repeat the detail here. The harmonics are 
essentially suppressed for Vq = 0.003, with a very nar- 
row dip at X = 1.2. But at Vq = 0.012, the harmonics 
are barely emerging, and there is a new dip structure 
developing at X = 1.4. At Vq = 0.024, the harmonic am- 
plitudes become very large. In the short potential range 
case, as shown in Fig. 3(b), the dip structure is narrower 
for the lower Vo values. There is no harmonic structures 
and the dip location shift only slightly towards the lower 
X region in the larger Vo case. Even though at X = 1.2, 
we have a ~ 0.19A, the result does not necessary reduce 
to the delta-profile result. We point out that only the 
Vo = 0.02 case in Fig. 3(b) corresponds to a delta-profile 
result for Vd = Vq(1 = 0.2. To demonstrate this , we 
show in Fig. 3(c) the a — 10, and Vo = 0.036 result 
(solid curve), and compare with the delta-profile result 
for Vd = 0.36 (the dotted curve). The two curves are very 
different. The most important different is that the dip lo- 
cation in the delta- profile case is shifted to X « 1.1 while 
the dip remains essentially unchanged for the FRTM po- 
tential. 

To further establish the conditions when the FRTM 
potential result reduces to a delta-profile result, we show 
in Fig.[||the FRTM potential result for a small V a = 0.2 
value, while varying a. We recall that at X = 1.2, A ~ 53. 
The a = 10 result (solid curve) is identical with the delta- 
profile result. However, for a — 20, and 30, the FRTM 
results are different from that of a delta-profile. As a 
is increasing, the dip at X — 1.2 becomes narrower and 
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FIG. 6: G as a function of X for V = 0.012, and u = 0.014. 
The interaction ranges are: "a" a — 150, "b" a — 100, "c" 
a — 50. The curves are vertically offset for clarity. 
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FIG. 7: G as a function of X for a = 200, and Vo = 0.012. 
The frequencies are: "a" lu = 0.014, "b" u = 0.028, "c" 
uj = 0.042. The curves are vertically offset for clarity. 



steeper, while the peak structure at X — 1.4 becomes a 
very shallow dip. This shows clearly that a <C A must 
also be satisfied. Thus the correct conditions for a FRTM 
potential result to reduce to that of a delta-profile result 
are a<A, and Vba <C 1. 

In Figs. 5 and 6, we fix Vo at 0.02, and varying a. Fig. 
5 is for the short potential range, with a < A. Within this 
short a regime, the dip at X = 1.2 becomes wider and 
shallower, as a is increasing. The dip at X = 1.4, how- 
ever, becomes more pronounced as a is increasing. These 
features indicate that high order processes such as the 
2uj processes have become significant. An electron that 
is trapped can then be excited back out of the QBS. Thus 
we find the general trend that when the 2oj processes are 
significant, the first dip will become shallower. Fig. 6, 
on the other hand, is for the long potential range cases. 
Harmonics appear in the large a regime. We note that 
the a = 50 result, which corresponds to the case a ~ A, 
exhibits the emerging effects of the harmonics, and the 
dip structures become very shallow. For larger a, the dip 
structures become larger again. 

In Fig. [7] we present the case for a = 200, Vo = 0.012, 
and varying ui. In curve a, the dip structures are sub- 
jected to the effect of the harmonics, since at the dip 
location, the harmonic amplitude is not that small. How- 
ever, for curves b, and c, the harmonic amplitudes are 
very small near the location of the dip. We point out also 
that the electron wavelengths A near the dip structure 



in the curves b, and c, are 37.6, and 30.7, respectively. 
Thus, for example, in curve c, a ~ 6.5A, and we are in 
the very long range regime. The QBS features are still 
very clear. 

In Fig. [SI we plot the time-average of the spatial prob- 
ability density for a = 200, V = 0.012, uj = 0.028, and 
X = 1.399. The dip location is at X = 1.4. Our choice 
of the parameters is near that for the occurrence of the 
QBS. The probability density shows the evanescent na- 
ture of the trapped electron. Furthermore, the higher 
probability density near the two edges of the FRTM po- 
tential shows that the QBS processes take place more 
frequently within a distance of order A from the edges 
of the profile, at which spatial variation in the potential 
occurs. 

Finally, in Fig. EI w e present the one sideband approx- 
imation results, together with our nonperturbative re- 
sults. The deviation is large in the X < 1.2 region, which 
corresponds to the case when the evanescent modes play 
an important role. The deviation is larger for larger Vo, 
as expected. This shows that the one sideband approxi- 
mation is not appropriate and, in fact, the conservation 
of current condition is usually not satisfied. 
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FIG. 8: Time-averaged probability density (| ip(x,t) | 2 ) as 
a function of longitudinal position £. The parameters are 
a = 200, Vo = 0.012, and X = 1.399. 
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FIG. 9: G as a function of X for a = 10, and u = 0.014. Non- 
perturbative result (dashed curve) and one-sideband approx- 
imation result (solid curve), (a) Vo = 0.01, (b) Vq = 0.005. 



tic and coherent. We find QBS features in all poten- 
tial ranges, including the long as well as the short range 
regime. The dip structures associated with the QBS oc- 
cur when \i is at mfko above the threshold of a subband 
edge. We find also that the inelastic processes occur more 
likely in the region when the potential profile varies spa- 
tially. In addition, from our results which we have not 
shown here, we find that a one-sideband approximation, 
in general, violates the conservation of current require- 
ment. 

We have presented arguments for the implications of 
our abrupt-profile FRTM potential results to that of a 
smooth-profile FRTM potential. We summarize and sup- 
plement our arguments, in the following, in light of the 
similarity and the difference between the two potential 
profiles. The abrupt- profile FRTM potential is similar 
to the smooth-profile FRTM potential in that they both 
break the longitudinal translational invariance. This al- 
lows the electrons to absorb or emit energy in units of 
fru, for arbitrary u>. Consequently the electrons can make 
transitions, via inelastic processes, to the QBS just be- 
neath a subband edge, giving rise in G to dip structures. 
A conclusion from this similarity between the profiles is 
that inelastic processes leading to QBS features are per- 
mitted in both of the profiles. 

These two potential profiles are different, however, in 
that the abrupt-profile potential introduces additional 
multiple scatterings between the two abrupt edges of the 
potential, and gives rise to harmonics in G. Due to these 
additional multiple scatterings, the electrons effectively 
stay longer within the abrupt-profile FRTM potential re- 
gion than within the smooth-profile FRTM potential re- 
gion. As a result, the QBS features in the former poten- 
tial profile might be perturbed, either being enhanced or 
suppressed. By proposing a small-harmonic-magnitude 
criteria, we attempt to look at cases where the perturba- 
tion from the harmonics to the QBS features is expected 
to be small. In these cases, the features of the two profiles 
are expected to be similar qualitatively. 

In conclusion, the CIS and the QBS features are found 
in a NC, acted upon by an abrupt-profile FRTM poten- 
tial. The features are argued to exhibit in the case of 
a smooth-profile FRTM potential, in particular, and is 
expected to affect the time-dependent properties of NC, 
in general. Further study, however, is needed to attain a 
better understanding. 



IV. CONCLUSION 
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